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Abstract We present a sample of 30 well observed Post Common Envelope Binaries (PCEBs). Deriving the cooling 
age of the white dwarfs, we show that the PCEB population is dominated by young systems. Having calculated 
the orbital evolution of the systems under the assumption of two different prescriptions for the angular momentum 
loss, we find that most of the systems have not yet completed a significant fraction of their PCEB life time. We 
therefore predict the existence of a large population of old PCEBs containing cold white dwarfs (T c s <15 000 K). 
Our calculations show that nearly half of the PCEBs in our sample will evolve into a semi-detached configuration 
and start mass transfer in less than the Hubble-time. These systems are thus representative for progenitors of 
the current CV population. Only one of them ( V471 Tau) will evolve into a long-period (P or b iJ 4 h) CV, and a 
(f) , significant fraction of the systems will start mass transfer in the period gap. Having estimated the distances of 

■ the PCEBs in the sample, we derive a space density of Ppceb ~ 6 — 30 x 10 -6 pc -3 , depending on the assumed 
angular momentum loss prescription. Taking into account the evolutionary time scales we compute a lower limit 
for the CV space density, predicted by the currently known PCEB population of pcv ii 10~ 5 pc -3 . Finally, we 
discuss possible observational selection effects and conclude that the observed PCEB population is probably highly 
incomplete. 

|' Key words, accretion, accretion discs - binaries: close - novae, cataclysmic variables. 
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c/3 , 1. Introduction magnetic braking (Rappaport et al. 1983; Paczynski & 

d ■ Sienkiewicz 1983; Spruit & Ritter 1983; see King 1988 for 

- In cataclysmic variables (CVs) (Warner, 1995, for an en- a review) _ In briefj the concept of this theory is that the 

M-j . cyclopaedic review) a white dwarf accretes material from evolution of CVs is divided into two main phaseSj depe nd- 

■ its Roche-lobe filling secondary, typically a K- or M-dwarf. ing on thc prev ailing angular momentum loss mechanism. 

The standard picture for the formation of CVs assumes Stellar magnetic braking dominates in CVs whose Roche- 
that the progenitor systems were moderately wide binaries l be filling donor stars still have a radiative core, which 
consisting of an intermediate mass star with a low mass i s the case for orbital periods P orb > 3 h. Once that the 
companion. Once the more massive star evolves to giant donor s t ars become fully convective at P orb ~ 3 h, mag- 
dimensions and fills its Roche-lobe, runaway mass trans- net i c braking ceases. For P orb < 3 h gravitational radiation 
fcr onto the less massive star starts and the systems en- takes over as a much less efficient angular momentum loss 
ters into a Common Envelope phase (CE). Friction within mechanism, resulting in longer evolution time scales. As 
the envelope extracts angular momentum, which tightens a consequence of the high mass loss rate in the magnetic 
the orbit. Post common envelope binaries (PCEB) with braking regime, the donor stars in CVs with P orb > 3h are 
orbital periods of a few days or less may evolve into a somewhat expanded. The mass loss rate decreases when 
semi-detached CV configuration through orbital angular magnetic braking ceases at P orb ~ 3 h, and the donor star 
momentum loss. Plausible angular momentum loss agen- reacts by relaxing to its thermal equilibrium configuration 
cies are gravitational radiation and - much more efficient - w ith a radius that is smaller than its Roche-lobe radius, 
magnetic braking. As a result, the mass transfer shuts off completely, and the 

The standard paradigm describing the evolution of CV becomes an inactive detached white dwarf /dM binary 

CVs after the onset of mass transfer is known as disrupted that evolves towards shorter periods through emission of 

gravitational radiation. At P orb ~ 2 h, the secondary fills 

Send offprint requests to: M.R. Schreiber, e-mail: again its Roche volume and re-starts the mass transfer 
mschrei@astro . u-strasbg . f r 
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(though at a much lower rate than in the long-period 
CVs). 

The main merit of the disrupted magnetic braking 
model is that it can successfully explain the period gap, 
i.e. the statistically significant paucity of known CVs with 
orbital periods in the range 2 — 3h. However, a number of 
predictions of this standard model are in strong disagree- 
ment with the observations: (1) the predicted minimum 
orbital period is 65 — 70min, 10% shorter than the ob- 
served value (Kolb & Baraffe, 1999); (2) CVs should spend 
most of their lifetime near the minimum period, increas- 
ing the discovery probability for such systems. Population 
syntheses predict a significant accumulation of systems 
near the minimum period, which is not observed (Kolb 
& Baraffe, 1999); (3) whereas all population syntheses 
predict that <~ 95% of the entire CV population should 
have orbital periods < 2h (e.g. Howell et al., 1997), sim- 
ilar numbers of CVs are observed above and below the 
period gap; (4) while the population models predict a 
space density of a few 10~ 5 to a few 10~ 4 pc~ 3 (de Kool, 
1992; Politano, 1996), the space density derived from 
the currently known sample of CVs is only several few 
1(T 6 pc~ 3 (e.g. Downcs, 1986; Ringwald, 1996). If the pop- 
ulation models are correct, we have identified so far only 
a small fraction (~ 1 — 10%) of the existing CV popu- 
lation (Gansicke et al., 2002); and finally (5) there is no 
observational evidence for a discontinuous change in the 
spin-down rate due to magnetic braking between late- type 
field stars that are fully convective and those that have a 
radiative core. 

It is apparent that detailed populations studies of the 
CV progenitors (i.e. PCEBs) are extremely important for 
a global understanding of CV formation and evolution, 
and extensive theoretical analyses of these systems have 
been performed in the past (Ritter, 1986; de Kool, 1992; 
de Kool & Ritter, 1993; King et al., 1994; Politano, 1996). 
However, the relatively small number of known PCEBs 
limited so far the comparison of these studies with obser- 
vations. Throughout the last decade, a number of addi- 
tional PCEBs have been discovered, significantly improv- 
ing the statistical properties of the known PCEB popula- 
tion. 

In this paper, we analyse the properties of a sample of 
well-observed PCEBs and discuss possible implications for 
the PCEB/CV evolution. In Sect. 2, we briefly summarise 
recent alternatives/additions to the standard CV evolu- 
tion theory outlined in the Introduction. The different an- 
gular momentum loss prescription that have been used in 
the context of CV evolution are described in Sect. 3. We 
introduce our sample of PCEBs in Sect. 4, and discuss the 
past and future evolution of these stars in Sect. 5 and 6, re- 
spectively. Section 7 provides the distances to the PCEBs 
in our sample. In Sect. 8, we compute the space density of 
PCEBs, and use this result to estimate the space density 
of the present-day CV population. The effects of obser- 
vational biases are discussed in Sect. 9. Finally, Sect. 10 
summarises our findings. 



When discussing the properties of the present-day 
PCEB population in the context of the present-day CV 
population, specifically in Sect. 6 and 8, we will use the 
term pre- cataclysmic variables (pre-CVs) to denote those 
systems which can be regarded as representative for the 
progenitors of the current CV population - i.e. PCEBs 
that evolve into a semi-detached configuration in less than 
the Hubble-time 1 . 



2. Recent modifications of the "standard 
scenario" of CV formation and evolution 

The "standard scenario" of CV formation and evolution 
outlined in the introduction has essentially remained un- 
changed during the last two decades, even though being 
challenged with a number of alternative suggestions (Livio 
& Pringle, 1994; King & Kolb, 1995; Clemens et al., 1998; 
Kolb et al., 1998). Very recently, two far-reaching modifi- 
cations for the standard scenario have been proposed. 

Sills et al. (2000) and Pinsonneault et al. (2002) pre- 
sented theoretical models of the angular momentum evo- 
lution of low mass stars (0.1 — 0.5 M ) and compared 
their models to rotational data from open clusters of dif- 
ferent ages to infer the rotational history of low-mass 
stars and the dependence of initial conditions and rota- 
tional evolution on mass. The studies of Sills et al. (2000), 
Pinsonneault et al. (2002), and Andronov et al. (2003) 
have two important consequences for the theory of CV 
evolution. On one hand, angular momentum loss via mag- 
netic braking is less efficient above the fully convective 
boundary than in the standard CV model, and as a re- 
sult the evolution time scale of PCEBs containing a sec- 
ondary with a radiative core is expected to be 2 orders 
of magnitude longer than in the standard model. On the 
other hand, the observed angular momentum loss proper- 
ties show no evidence for a change in behaviour at the fully 
convective boundary. Magnetic braking remains, hence, an 
important angular momentum loss mechanism in PCEBs 
with a fully convective secondary, and these systems evolve 
faster than in the standard scenario. Due to the increased 
angular momentum loss below the gap, the theoretical 
orbital minimum increases, which is in better agreement 
with the observations (Kolb & Baraffe, 1999; Patterson, 
1998; King et al., 2002). 

However, the evolution time scale of CVs above the 
gap is much longer than in the standard scenario, and 
the existence of the observed period gap is tentatively ex- 
plained by Pinsonneault et al. (2002) by two separated 
populations of CVs with differently evolved donor stars. 
One problem with this modification of the standard CV 
evolution scenario is that the reduced angular momentum 
loss above the period gap predicts mass transfer rates that 

1 Obviously also PCEBs with evolutionary time scales longer 
than the Hubble-time will form CVs, but they differ from the 
progenitors of present-day CVs. This difference is the only rea- 
son for our definition of "pre-CVs" . 
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are significantly lower than the values derived from the ob- 
servations. This issue is discussed in more detail in Sect. 3. 

King & Schenkcr (2002) and Schcnkcr & King (2002) 
postulate that the PCEB evolution into a semi-detached 
configuration takes much longer than in the standard sce- 
nario. Without modifying the standard angular momen- 
tum loss prescription, such a situation would arise if the 
frictional angular momentum loss during the CE is less 
efficient, the PCEBs exit the CE with wider orbits, and 
evolve into the semi-detached CV configuration only on 
longer time scales than typically assumed. If the time 
scale to initiate mass transfer is of the order of the galac- 
tic age (~ I0 10 yrs), both the minimum orbital period 
problem and the space density problem are solved as a 
much smaller number of CVs has been produced to date 
compared to the standard scenario, and the present-day 
CV population has not yet reached the minimum period. 
Finally King & Schenker (2002) speculate that the longer 
evolution time scale for PCEBs proposed in their model 
might be related to a reduced magnetic braking as sug- 
gested by Andronov et al. (2003). However, if the assumed 
birth rate of the PCEB progenitor binaries remains un- 
changed, a longer PCEB phase must result in a very large 
population of old PCEBs. In the following sections, we 
analyse the properties of the known PCEB sample and 
discuss the results in the framework of both the standard 
CV evolution and the context of a revised magnetic brak- 
ing prescription. 

3. Angular momentum loss 

Angular momentum loss (AML) drives the evolution of 
the binaries and hence, our understanding of AML is the 
main ingredient of CV and PCEB evolution. It is believed 
that PCEBs as well as CVs lose angular momentum via 
gravitational radiation and magnetic braking. 

The rate of AML due to radiation of gravitational 
waves from Einstein's quadrupolc formula is 



Jgr = - 



32 G 7 / 3 Mi 2 M scc 2 / 2tt 
5c 5 (Mi + M sec ) 2 / 3 \P~^ 



7/3 



(1) 



with Mi and M sec the mass of the primary and the sec- 
ondary. 

Unfortunately, the efficiency of magnetic braking 
is rather uncertain (see the Introduction and Sect. 2) 
and it may have been overestimated in the past (see 
Pinsonneault et al. , 2002) . In the following analysis of the 
currently known PCEB sample we use, therefore, two dif- 
ferent prescriptions of AML to calculate the evolution of 
our sample of PCEBs. 

(1) Disrupted magnetic braking until the secondary be- 
comes fully convective and magnetic braking ceases (fol- 
lowing the standard scenario of Verbunt & Zwaan, 1981; 
Rappaport et al., 1983). In this case, the AML is given by 



j vz = -3.8 x IO- 3O M sec i? 4 



Rs, 



R(7) 



2tt 



orb 



(2) 



The original version of the magnetic braking law corre- 
sponds to 7 = 4.0 (Verbunt & Zwaan, 1981). In the con- 
text of the CV evolution 7 = 2.0 is frequently used, and 
hence we adopt this value throughout this work. 

(2) The empirical AML prescription derived from open 
cluster data of single stars is (Sills et al., 2000), 



Js 



YV 



PT 



RsccMq 

M scc Rq 

-^{m^Rq- 



0.5 



0.5 



forw < w crit 

forw > Wcritj 



(3) 



where w C rit depends on the mass of the secondary. For 
stars with masses M > M Q a Rossby scaling matches the 
observations 



^crit — W cr ;tQ 



r(M scc ) ' 



(4) 



with t being the convective overturn time (Krishnamurthi 
et al., 1997). For lower masses the dependence of w cr i t on 
mass is found to be stronger than given by Eq. (4) (Sills 
et al., 2000). Throughout this work we adopted the values 
for o> cr it given by Andronov et al. (2003). Notice, for the 
systems in our sample we find u = 10 — f0 2 x w cr i t . 

In the following sections we distinguish between the 
two currently discussed AML prescriptions using the no- 
tation Jcmb and Jrmb when referring to the classical re- 
spectively the reduced AML prescription: 



Jcmb 
and 



Jvz 



Jgr 
Jgr 



for 
for 



M sec > 0.3M F 
M sen < 0.3M 



o 







Jrmb = Jspt + Jgr- 



(5) 



(6) 



Fig. 1 shows the AML as a function of P rb for PCEBs. 
Fig. 2 compares the angular momentum loss agencies 
Jcmb, Jgr, and Jrmb for CVs (see also Pinsonneault 
et al., 2002; Andronov et al., 2003). 

It is worth noting that the reduced magnetic braking 
prescription is not only subject to observational uncertain- 
ties (as stated by Andronov et al. (2003) themselves) but 
also unable to explain the mass accretion rates derived 
from observations for CVs above the orbital period gap 
(e.g. Patterson, 1984). This is a firm conclusion unless (I ) 
all nova-like systems represent short-lived high accretion 
states or (2) there exists an additional AML mechanism 
in these systems. Referring to (I) we note, the possibil- 
ity of irradiation induced mass transfer cycles (e.g. King 
et al., 1995) and that nova eruptions widen the mass trans- 
fer spectrum (Kolb, 2002). Concerning (2) we note that 
circumbinary disks in CVs have been suggested as rather 
efficient additional AML agents (Spruit & Taam, 2001; 
Dubus et al., 2002). Thus, while the reduced AML pre- 
scription has a problem to explain the "observed" mass 
transfer rates of CVs above the gap, it is currently not 
possible to rule out its validity because of the additional 
complications just mentioned 2 . 

2 In this context we note that it might be easier to test our 
understanding of AML using a large (not yet identified) sample 
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Figure 1. Angular momentum loss in PCEBs as a func- 
tion of the orbital period. We assumed M scc — O.4M and 
Mi = 0.6 Mq. Apparently the revised magnetic braking pre- 
scriptions is of similar efficiency as the classical prescription for 
long orbital periods whereas it is only a few times the angular 
momentum loss due to gravitational radiation for short orbital 
periods. 




Figure 2. Angular momentum loss in CVs. As in Fig. 1 
we used Mi — O.6M (see also Pinsonneault et al., 2002; 
Andronov et al., 2003). The dotted lines correspond to Eq. (5) 
with two values for 7 (see Eq. (2)). The abrupt decrease dis- 
plays the shutoff of magnetic braking for M scc < O.3M . The 
dashed line represents angular momentum loss by gravitational 
radiation only whereas the solid line refers to the revised an- 
gular momentum loss prescription (Eq. (6)). 



Considering the described uncertainties of the AML in 
close binaries, we discuss the current PCEB population in 
the context of both prescriptions, Jcmb and Jrmb- 
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Figure 3. The ZAMS mass radius relation given by Politano 
(1996) and the positions of the PCEB secondaries for which 
an observationally determined radius is available. The arrows 
indicate that the radii of the secondaries in these planetary 
nebulae are larger due to recent accretion onto the secondary 
during a CE phase. Before becoming a CV the secondary will 
have contracted to its main-sequence radius. 

4. PCEBs: The present 

We have compiled for our present study a sample of 
PCEBs with well-established characteristics from various 
literature sources, including the catalogues of Ritter & 
Kolb (1998) and Kubc et al. (2002). We exclude from our 
sample systems with M scc > Mi (such as V651 Mon), sys- 
tems with P or b ^ 2d (such as Feige 24), and systems 
with sub-giant secondary stars (such as FF Aqr) . Table 1 
lists the properties of our PCEB sample, which consists 
of 18 systems containing a white dwarf primary and 12 
systems containing an sdOB primary (of which four are 
the central stars of planetary nebulae). These 30 systems 
represent the currently known population of well-observed 
PCEBs which will evolve eventually into a semi-detached 
configuration and turn into a CV and for which reliable 
observational data exists. 

It is interesting to note that the radii of the secondaries 
derived from observations (Table 1) are generally in good 
agreement with the ZAMS mass-radius relation given by 
Politano (1996). Exceptions are the secondary stars in the 
binary planetary nebulae BEUMa, KV Vel, V477Lyr, and 
UUSge whose radii exceed the ZAMS value by a factor 
of ^ 2 (Fig. 3). This effect can be explained by accre- 
tion of matter onto the secondary at a relatively high rate 
during the recent CE-phase. Numerical calculations of ac- 
creting low-mass stars have been performed by Prialnik & 
Livio (1985) for a fully convective star of M = 0.3M Q 
and Fujimoto & Iben (1989) for a low mass main se- 
quence star of M = 0.75M©. The calculations show that 
rather recent accretion at a rate of ~ lO _5 M yrs _1 to 

of PCEBs, as all the complications related to mass transfer are 
not present in these binaries. 
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Table 1. Sample of PCEBs with orbital periods P<2days, Mi > M scc and a main-sequence secondary star. 



Ubject 


alt. name 


bpl 


-Torb [dj 


Mi [Mq\ 


71 /T r 71 jT 1 
M S ec [Mq\ 


r> r D 1 

-Ksec [-K©J 


Jeff [K.J 


Kel. 


RRCae 


LFT 349 


WD 


0.304 


0.467 


0.095 


0.189 


7000 


1,2 


EGUMa 


Case 1 


WD 


0.668 


0.64 + 0.03 


0.42 ± 0.04 


0.40 % 


13125 ± 125 


3 


EC 13471-1258 




WD 


0.151 


0.77 + 0.04 


0.58 + 0.05 


0.42 


14085 ± 100 


4 


BPM 71214 




WD 


0.201 


0.77 + 0.06 


0.4 


0.38 % 


17200 ± 1000 


4 


HRCam 


GD448 


WD 


0.103 


0.41+0.01 


0.096 ± 0.004 


0.125 + 0.020 


19000 


5,6 


UZSex 


PG 1026+002 


WD 


0.597 


0.65 + 0.23 


0.22 + 0.05 


0.25 + 0.03 


17600 ± 2000 
19900 ± 330 


7 

8 


BPM 6502 


LTT 3943 


WD 


0.337 


0.5 + 0.05 


0.16 + 0.09 


0.204 % 


20640 ± 200 
21380 ± 258 
20311 + 532 


9 
2 
10 


HZ 9 




WD 


0.564 


0.51+0.1 


0.28 ± 0.04 


0.296 % 


20000 ± 2000 
17400 


11,12,13 
this work 


MS Peg 


GD245 


WD 


0.174 


0.48 + 0.02 


0.22 + 0.02 


0.27 + 0.02 


22170 


14 


CCCet 


PG 0308+096 


WD 


0.284 


0.39 + 0.1 


0.18 + 0.05 


0.21 + 0.02 


26200 ± 2000 


7 


HWVir 


BD-07°4377 


sdB 


0.117 


0.48 + 0.09 


0.14 + 0.02 


0.180 + 0.011 


28488 ± 208 


15 


HS 0705+6700 




sdB 


0.096 


0.483 


0.134 


0.186 


28800 


16 


LM Com 


Ton617 f 


WD 


0.259 


0.45 + 0.05 


0.28 ± 0.05 


0.296 % 


29300 


17 


PC 1017-086 




sdB 


0.073 


0.5 


0.078±Ho6 


0.085 ± 0.04 


30300 


18 


V471 Tau 


BD+16°516 


WD 


0.521 


0.84 + 0.05 


0.93 ± 0.07 


0.98 + 0.10 


34500 ± 1000 


19 


NY Vir 


PG 1336-018 


sdB 


0.101 


0.5 


0.15 


0.205 + 0.01 


33000 ± 3000 


20 


AA Dor 


LB 3459 


sdO 


0.262 


0.330 ± 0.003 


0.066 + 0.001 


0.10 + 0.01 1 


42000 ± 1000 


21 


RE 2013+400 




WD 


0.706 


0.56 + 0.03 


0.18 + 0.04 


0.221 % 


49000 ± 700 
47800 ± 2400 


22 
23 


GKVir 


PG 1413+015 


WD 


0.344 


0.51+0.04 


0.1 


0.15 


48800 ± 1200 


24,25 


MT Ser 


PN A66 41* 


sdO,PN 


0.113 


0.6 


0.2 + 0.1 


0.2 + 0.1 


50000 


26,27 


IN CMa 


REJ0720-318 


WD 


1.26 


0.57 + 0.02 


0.39 ± 0.07 


0.375 % 


53000 ± 1100 


22 


NNSer 


PG 1550+131 


WD 


0.130 


0.57 + 0.04 


0.12 ± 0.03 


0.166 % 


55000 ± 8000 


28 


TW Crv 


EC 11575-1845 


sdO 


0.328 








95000 ± 40000 


29 


RE 1016-053 




WD 


0.789 


0.6 + 0.02 


0.15 + 0.02 


0.195 % 


55000 ± 1000 
56400 ± 1200 


22 
23 


UU Sge 


PN A66 63* 


sdO,PN 


0.465 


0.63 + 0.06 


0.29 ± 0.03 


0.53 + 0.02 


87000 ± 13000 
57000 ± 8000 


30,31,32 
32 


V477Lyr 


PN A66 46* 


sdOB,PN 


0.472 


0.51+0.07 


0.15 + 0.02 


0.46 ± 0.03 


60000 


33 


PN A66 65 


Abell 65 


sd?,PN 


1.00 








80000 


34,35 


KVVel 


LS 2018 


sdO,PN 


0.357 


0.63 + 0.03 


0.25 + 0.06 


0.402 ± 0.005 


77000 ± 3000 
90000 


36 
37 


HS 1136+6646 




WD 


1 








100000 


38 


BEUMa 


PG 1155+492 


sdO,PN 


2.291 


0.7 + 0.07 


0.36 ± 0.07 


0.72 + 0.05 


105000 ± 5000 


39,40 



Other frequently used designations: * Abell 41, 46 & 63 (to be used with care as the object numbers in Abell 1955 and Abell 
1966 are not consistent); f PG 1224+309. 

! The radius of the secondary is estimated using the mass-radius relation shown by Kudritzki et al. (1982) (their Fig. 8). 
% Radii of the secondary calculated by using the mass-radius relation given by Politano (1996). 

References: (1) Bruch 1999, (2) Bragaglia et al. 1995, (3) Bleach et al. 2000, (4) Kawka et al. 2002, (5) Marsh & Duck 1996, 
(6)Maxted et al. 1998, (7) Saffer et al. 1993, (8) Kepler & Nelan 1993, (9) Kawka et al. 2000, (10) Koester et al. 1979, 
(ll)Stauffer 1987, (12)Lanning & Pesch 1981,(13) Guinan & Sion 1984 (14) Schmidt et al. 1995, (15) Wood & Saffer 1999, 
(16)Drechsel et al. 2001, (17) Orosz et al. 1999, (18) Maxted et al. 2002, (19) O'Brien et al. 2001, (20) Kilkenny et al. 1998, 
(21) Rauch 2000, (22) Vennes et al. 1999, (23) Bergeron et al. 1994, (24) Fulbright et al. 1993, (25) Green et al. 1978, (26) Green 
et al. 1984, (27) Grauer & Bond 1983, (28) Catalan et al. 1994, (29) Chen et al. 1995 (30) Miller et al. 1976, (31) Pollacco & Bell 
1993, (32) Bell et al. 1994, (33) Pollacco & Bell 1994, (34) Bond & Livio 1990, (35) Walsh & Walton 1996, (36) Hilditch et al. 
1996, (37) Herrero et al. 1990, (38) Sing et al. 2001, (39) Wood et al. 1995, (40) Ferguson et al. 1999. 



~ lO _7 M yrs _1 can explain the observed expansion of 
the secondary. As the time scale on which the secondary 
contracts is shorter than the binary evolution time scale 
we use radii calculated from the mass-radius relation of 
Politano (1996) for the four planetary nebulae in our cal- 
culations (Sect. 5, and 6). 



A peculiar system is the Hyades PCEB V471 Tau, lo- 
cated in the upper right corner of Fig. 3. Apparently, the 
secondary is expanded with respect to the ZAMS mass ra- 
dius relation. Using radiometric measurements and ground 
based Doppler imaging O'Brien ct al. (2001) explain the 
oversized secondary in V471 Tau as the result of a large 
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fraction of the secondaries surface being covered by star 
spots, lowering its average effective temperature. In order 
to maintain the luminosity of a 0.93 Mq dwarf, the ra- 
dius of the star has, hence, to exceed the ZAMS value. It 
is worth to note that also the white dwarf in V471Tau 
represents an evolutionary paradox as it is the hottest 
(and therefore youngest) as well as the most massive white 
dwarf in the Hyades cluster. O'Brien et al. (2001) dis- 
cuss the possibility that V471 Tau is descended from a 
blue straggler, which would imply that its progenitor was 
initially a triple system. Nevertheless, despite its curious 
nature, V471 Tau is the best-observed PCEB, and, as de- 
scribed in detail in Sect. 6, it will be the PCEB from our 
sample that starts mass transfer at the longest orbital pe- 
riod. 

Final notes concern MTSer, UUSge, EC 13471-1258, 
BPM 71214, PC 1017-086 and AA Dor. The orbital period 
of MTSer might be as twice as long as given in Table 1, 
and the system may consist of two hot subdwarfs rather 
than a subdwarf plus cool companion (Bruch et al., 2001). 
The large discrepancy between the two values for the tem- 
peratures of the primary in UU Sge given in Table 1 results 
from the fact that Bell et al. (1994) could fit the observed 
light curve equally well by assuming either limb-darkening 
or limb-brightening. The given values may therefore be 
interpreted as upper and lower limits. The binary param- 
eter for EC 13471-1258 and BPM 71214 given in Table 1 
indicate these systems are extremely close to the semi- 
detached state. Therefore Kawka et al. (2002) suggested 
that at least EC 13471-1258 might also be a hibernat- 
ing nova instead of a PCEB. We additionally note that 
the radius obtained for the secondary in EC 13471-1258 
is essentially lower than predicted by the main-sequence 
mass radius relation which appears implausible and in- 
dicates rather large uncertainties in the values for M sec 
and/or R scc . Finally, the classification of the secondaries 
in PC 1017-086 and AADor is not unambiguous. It has 
been noted that in both systems the secondary may also be 
a brown dwarf rather than a low mass M-dwarf (Maxted 
et al., 2002; Rauch, 2000). 

5. PCEBs: The past 

5.1. The white dwarf cooling age 

The cooling of single white dwarfs has been modelled to 
a very high degree of precision allowing to derive ac- 
curate age estimates for field white dwarfs from their 
(observationally determined) effective temperatures and 
masses. Some uncertainties in the theoretical models re- 
main at the low-temperature end of the white dwarf lumi- 
nosity function, and are relevant primarily for the cosmo- 
chronological interpretation of the coldest (halo) white 
dwarfs (see the reviews by Fontaine et al. 2001; Koester 
2002). 

We determined the cooling ages for the PCEBs in our 
sample which contain a white dwarf primary by interpo- 
lating Wood's (1995) CO core evolution tracks for the ef- 



fective temperatures and masses listed in Table l 3 . The 
cooling age estimates are given in Table 2, and are shown 
in Fig. 4 along with Wood's cooling tracks. The bottom 
panel of Fig. 4 displays the age distribution of our PCEB 
sample. It is evident that the currently known population 
of PCEBs is strongly dominated by systems younger than 
5 x 10 8 yrs, with the exception of RR Cae, which has a cool- 
ing age of t C ooi ~ 10 9 yrs. The mean cooling age of the 18 
systems in the sample containing a white dwarf primary 
is Icooi ~ 1-4 x 10 s yrs. Considering that the sdOB, sdO, 
and PN systems in our sample arc white dwarf progen- 
itors, the mean age of the entire PCEB sample is even 
lower. The accuracy of our cooling age estimates depends 
(a) on the uncertainties in the white dwarf masses and (b) 
on the core composition of the white dwarfs. We discuss 
both issues in more detail below. 

(a) The white dwarf masses in PCEBs can be measured 
to a rather high level of precision, as several independent 
methods are available in these systems (fitting model spec- 
tra to the Balmer lines, dynamical measurements from the 
radial velocity variations of the white dwarf and its com- 
panion, measurement of the gravitational redshift of pho- 
tosphcric lines from the white dwarf). The errors on the 
white dwarf masses extracted from the literature (Table 1) 
are completely negligible in the context of our cooling age 
estimates. Even if we assume that the published errors 
are underestimated, and that the true errors in the white 
dwarf masses are of the order ±0.1 — O.2M , the corre- 
sponding uncertainties in the derived cooling ages are of 
the order ~ 2, which is still irrelevant for the comparison 
to the binary evolution time scales. 

(b) The evolution of a white dwarf progenitor in a bi- 
nary system may affect the core composition of the white 
dwarf. If the first Roche lobe overflow phase occurs prior 
to helium ignition in the primary star (de Kool & Ritter, 
1993; Iben & Tutukov, 1993) the resulting white dwarf 
will be of low mass and have a helium core. Consequently, 
evolution models for CVs predict that systems with a 
low mass primary (Mi < O.5M ) should contain a he- 
lium core white dwarf (de Kool & Ritter, 1993; de Kool, 
1992; Politano, 1996; Howell et al., 2001). Careful observa- 
tions of seven previously considered low mass single white 
dwarfs revealed that five of them are indeed close binary 
systems (Marsh et al., 1995). 

Comparing the cooling models for CO white dwarfs 
(Wood, 1995) and He white dwarfs (Althaus & Benvenuto, 
1997; Driebe et al., 1998, 1999) shows, however, that the 
evolution of CO and He white dwarfs differs noticeably 
only for very low-mass stars, where residual hydrogen shell 
burning dominates over the gravothermal energy release. 
Inspecting Figs. 2 and 3 of Driebe et al. (1998) shows 
that only for a single system from our sample (Table 1) 
the cooling age derived from Wood's (1995) CO tracks 

3 Note that the same method can not be applied in CVs, as 
accretion (re)heats the white dwarf. As a result, white dwarfs in 
CVs are typically hotter than field white dwarfs of comparable 
age and mass (Gansicke, 2000; Townsley & Bildsten, 2002). 
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Figure 4. Cooling tracks for non-accreting white dwarfs after Wood (1995) for different masses of the white dwarf. The position 
of the currently known PCEBs (Table 1) are obtained by interpolating between the cooling tracks for different masses of the 
primary. The bottom panel shows the age of the PCEB population in a histogram with logarithmic bins (shaded, N) and the 
relative cumulative distribution (solid line, E). 



may be significantly wrong: if the white dwarf in RRCae 
(0.467 M , 7000 K) contains a He core, then its cooling age 
is ~ 4 — 5 x 10 9 yrs, compared to 1.2 x 10 9 yrs for a CO 
core. We conclude that the effect of different core compo- 
sitions is negligible in the context of the present paper, 
but has to be taken into account in future analyses if ad- 
ditional PCEBs containing low mass (Mi SS 0.4 M©) and 
cool (T e ff Si 20000 K) white dwarfs are found. 

A final note concerns the possible effect of re-accretion 
of CE material, which may alter the composition of the 
white dwarf envelope in a PCEB with respect to that of 
a single white dwarf. Differences in envelope/atmosphere 
composition will affect the cooling of the white dwarf only 
at the low-temperature end of the luminosity function. 
While this effect is important for the cosmochronological 
interpretation of the coldest and oldest halo WDs (e.g. 
Koester, 2002), it is irrelevant for our PCEB sample. 

5.2. The orbital period at the end of the common 
envelope phase 

During the detached post CE-phase the masses of the com- 
panions remain essentially constant, i.e. M wc j = M scc = 0. 
With Kepler's third law this leads to the well known rela- 
tion 



J 
J 



P 

3P' 



(7) 



Based on the the cooling age of a given PCEB derived 
from cooling tracks (Fig. 4, Table 2) we can approximate 
the orbital period at the end of the CE phase by integrat- 
ing Eq. (7). We have done this calculation for both the 
standard scenario (Eq. (5)), and for the empirical AML 



prescription by Sills et al. (2000) (Eq. (6)). In the stan- 
dard scenario magnetic braking is assumed to be much 
stronger than gravitational radiation for M scc > 0.3M© 
and, hence, we can simply integrate Eq. (7) using either 
J = Jgr or J = J V z depending on the mass of the sec- 
ondary. For gravitational radiation only we get: 

5c 5 (Mr + M scc ) 1 /3 



^cool 



256G 5 / 3 (2tt) 8 / 3 



X (PcE ~ ^orb ^ )j 

(see also Ritter, 1986). For systems with M 
we assume J = Jyz an d derive: 

2.63 1 29 G 2 / 3 Mi " • 
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Considering the empirical angular momentum loss pre- 
scription of Sills et al. (2000) the situation is somewhat 
more complex as the two angular momentum loss agen- 
cies can be of the same order of magnitude (see Fig. 1) and 
we can not neglect gravitational radiation. Nevertheless, 
there exists also an analytical solution for the integral of 
Eq. (7) using j = j R uB- 



^cool 



3(1 + aP 4 / 3 ) 



4a(P- 5 / 3 + p-l/3)pl/3 
3(1 + aP 4 / 3 ) ln(l + aP 4 / 3 ) 
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Having derived the cooling age of the PCEBs we solve 
Eqn. (8), (9), and (10) to obtain the orbital period im- 
mediately after the CE phase (Pce)- Again, the re- 
sults are given in Table 2. Not surprisingly, P CE is very 
close to the present orbital period of the systems as we 
found most of the PCEBs to be quite young. For sys- 
tems with t coo i<10 7 yrs we get essentially P or b — Pce- 
Clearly, the assumed AML prescription makes a signifi- 
cant difference only for rather old PCEBs (i.e. RRCae, 
in our sample) or for systems not too young and with 
M scc > 0.3M© (EGUMa, V471Tau, EC 13471-1258, 
BPM 71214). In these cases we give in Table 2 both re- 
sults where the second line corresponds to J = Jrmb- 
Comparing the obtained distribution of orbital periods 
after the CE-phase (Pce) with theoretical predictions 
(Fig. 9 in de Kool, 1992), we find reasonable agreement: 
for our systems we get 0.4 > log(PcE) > —1.1 with <~ 60% 
in the range of —0.1 > log(PcE) > —0.6. 

5.3. The total age of the PCEBs 

The cooling age of the primary represents just the age 
of the systems after the end of the CE. To approximate 
the total age of the binary, i.e the time since it appeared 
on the main sequence we use simple analytical fits to the 
stellar evolution, Roche geometry, and Kepler's third law 
to obtain the nature of the initial main sequence binary. 
Combining Eqn. (35,36), and (43)-(49) of Politano (1996) 
with the standard prescription of the binary shrinkage 
during the CE-phase (Eq. (8) in de Kool, 1992) gives the 
mass of the initial primary and, hence, the nuclear evolu- 
tion time. This is, of course, just a rough approximation 
because an uncertain fraction «ce of the gravitational 
binding energy that is released when the secondary and 
the primary spiral together goes into the ejection of the 
envelope. However, we assume q:ce = 1 and A = 0.5 (see 
de Kool, 1992) which gives results in reasonable agree- 
ment with the initial-to-final mass relation for CO-white 
dwarfs (Weidemann, 2000). The estimates obtained for the 
masses of the progenitors of the primaries (Mi jP ce) and 
the resulting evolution time since the binary appeared on 
the main sequence (i eV oi) are also given in Table 2. For 
PCEBs containing a low mass white dwarf t cvo \ can be 
of the order of 10% of the Hubble time (see Table 2) and 
therefore is not generally negligible. 

6. PCEBs: The future 

6.1. How long before the onset of mass transfer? 

In the last section we showed that the vast majority of 
the known PCEBs are young objects. Here, we calculate 
the orbital periods that these binaries will have when they 
turn into CVs as well as the time scale for this evolution. 

Knowing Mi and M scc of a (detached) PCEB (Tabic 1) 
the orbital period of the corresponding semi-detached con- 




' log(Porb [d]) 

Figure 5. The orbital period distribution of the currently 
known CV population (from Kube et al. 2002). Grey-shaded 
is the orbital period gap. The vertical lines on the top indi- 
cate the periods at which the PCEBs from Table 1 will enter 
a semi-detached configuration and start mass transfer. Those 
PCEBs that will evolve into CVs within less than a Hubble 
time (Table 2) are considered to be genuine pre-CVs. The two 
PCEBs with P s d below the minimum orbital period of CVs at 
~ 80min are the two systems in our sample which possibly 
have a brown dwarf secondary AADor and PC 1017-086 (see 
Sect. 4). 

figuration, P s d follows from Roche geometry and Kepler's 
third law (Ritter, 1986): 

/ Rscc 3 \0.5 

Psd = UM S cc(l + M 1 /M S0C )(P L / a )3j ■ 

Here Pl denotes the volume radius of the Roche-lobe of 
the secondary and a the binary separation. We use the 
formula of Eggleton (1983) to relate these quantities to 
the mass ratio of the system. The resulting P s d for the 
PCEBs in our sample are listed in Table 2 and are com- 
pared to the orbital period distribution of the currently 
known CV population in Fig. 5. It is apparent that only a 
small number of systems (~ 22%) of the currently known 
PCEB population will evolve into CVs above the period 
gap (P s d > 3h), whereas ~ 30% will evolve into contact 
at P s d < 2 h. A somewhat surprising result is that almost 
half (~ 48%) of the currently known PCEBs will start 
their CV life in the period gap (2h< P s d < 3h). The 
small number of progenitors of long-period CVs will be 
discussed in Sect. 9. 

Knowing P s d we now calculate the time it will take a 
PCEB to become a CV (i s d), again by integrating Eq. (7) 
assuming either AML according to the "standard sce- 
nario" (j = J C mb) or Sill's et al. (2000) empiric AML 
prescription (J = Jrmb)- The solution for i s d can be de- 
rived by replacing i coo i, P orb , Pce with i sd , P sd , P or b in 
Eqn. (8-9). It is worth to note that although the PCEBs 
evolve on a time scale comparable or even longer than 
10 10 yrs, the expansion of the secondary due to its nuclear 
evolution is negligible for all the PCEBs in our sample 
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Figure 6. The PCEB evolution of HRCam and EG UMa as- 
suming two different prescriptions for the AML, Jcmb and 
Jrmb- As the mass of the secondary in HRCam is signifi- 
cantly below 0.3Mq the system evolves slower into contact if we 
assume AML according to the standard scenario. Conversely, 
EG UMa turns into a CV much faster in the standard sce- 
nario because of the more efficient magnetic braking for M scc > 
0.3M Q . 
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Figure 7. The cooling tracks of Wood (1995) and the loca- 
tions in the age/temperature plane where mass transfer will 
initiate for the PCEBs with a short PCEB lifetime (t CO oi + 
t sd < 6 x 10 9 yrs). The filled squares correspond to Jcmb the 
open squares to Jrmb- Notice, most of the PCEBs will stay 
more than 10 10 yrs in the PCEB state and, hence, will have 
T c ft < 6000 K when becoming semi-detached. 



and, hence, not taken into account in our calculations. 
Inspecting the resulting t s( j (Table 2), it becomes evident 
that most of the PCEBs have completed only a small frac- 
tion /pceb of their predicted PCEB-lifetime. On average 
the total evolution time for the PCEBs is ~ 3.4 x 10 10 yrs 
independent on the assumed AML prescription 4 . 

Assuming J = Jcmb, four systems (i.e. EG UMa, 
EC 13471-1258, BPM 71214, V471 Tau) have / PC eb > 
0.1. EC 13471-1258 is expected to be very close to the 
onset of mass transfer, /pceb > 0.995, and one may spec- 
ulate that this system is an old nova rather than a PCEB 
(Kawka et al., 2002, see also Sect. 4). 

As expected, i s d of an individual system depends 
strongly on the assumed AML prescription. Following 
Andronov et al. (2003) we obtain shorter evolution times 
(by a factor of 1.5 — 4) compared to the standard sce- 
nario for systems with M scc < 0.3M Q , whereas for sys- 
tems with Msec > 0.3 Mq i s d is significantly longer (by a 
factor of 3 — 30 depending on the present orbital period of 
the system). The differences in t s d, depending on the use 
of Jcmb or Jrmb are illustrated for the systems HR Cam 
and EG UMa in Fig. 6. 

Figure 7 shows that most of the PCEB white dwarfs 
will be very cool when the secondary is expected to fill its 
Roche-lobe and mass transfer starts. Only for JcMBwe get 
one system, (i.e. V471 Tau, the filled square in the upper 
left corner of Fig. 7) in which the mass transfer starts while 
the primary still can be considered a "hot" white dwarf. 

6.2. Representative progenitors of the current CV 
population 

Analysing PCEBs is not only interesting in itself, it is also 
important in the context of CV formation and evolution 
(see Sect. 1 and 2). Discussing the properties of our sam- 
ple of PCEBs in the framework of the present-day CV 
population requires to select those systems which are rep- 
resentative for the former progenitors of the current CV 
population. Only PCEBs which will evolve into CVs in 
less than the Hubble time (assuming t = 1.3 x 10 10 yrs 5 ) 
satisfy this condition. To distinguish these systems from 
the PCEBs which are different from the former progeni- 
tors of the present CVs we refer to them as pre- CVs (see 
also Sect. 1). 

Applying this selection criterion, we find 14 (16) pre- 
CV candidates in our sample (Table 2) when assuming 
Jcmb (Jrmb)- Although i s d depends strongly on the as- 
sumed AML for some individual systems (e.g. EG UMa, 

4 Notice, Jrmb is smaller than Jcmb only above the gap 
otherwise Jrmb is more efficient. 

5 Recently Ferreras et al. (2001) set new constraints on the 
age of the universe combining the results obtained by cos- 
mochromology, stellar population synthesis, and mapping of 
the peaks in the microwave background. They derived an age 
of 13.2±o;g Gyrs. Note also that Cayrel et al. (2001) obtained 
a stellar age of 12.5 ± 3 Gyrs from uranium decay. Throughout 
this work we therefore assume that star formation in the 
Galaxy began 13 Gyrs ago. 
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Table 2. The properties of the PCEBs calculated in Sect. 5-7. If there is no reference given in the Table, the distance is 
obtained from the references listed in Table 1. The time until the mass transfer is expected to start (t s a) is calculated using 
Eq. (8), (9) (CMB) and Eq. (10) (RMB). The total time until the system becomes semi-detached, t to t is given by the sum of 
the present cooling age, t C ooi, the estimated MS lifetime of the primary, t evo \, and t s d- PCEBs with t to t less than a Hubble time 
are marked in bold. For the systems RRCae, EG UMa, EC 13471-1258, BPM 71214 and V471 Tau we find the properties of 
the progenitor of the PCEB (orbital period Pce, primary mass Mi iP ce and, t e voi) depending on the assumed AML (see text). 
In these cases, the second line gives the results assuming reduced magnetic braking (RMB) whereas the first line corresponds 
to the classical magnetic braking assumption (CMB). /pceb gives the mean fractional PCEB life time that the system has 
already passed through. The orbital period at which the PCEB starts mass transfer P sd is independent on the assumed angular 
momentum loss. 
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(1971), (6) Maciel (1984), (7) Cahn et al. (1992), (8) Vennes & Thorstensen (1996), (9) Wood & Marsh (1991), (10) Walton 
ct al. (1993), (11) Pottasch (1996). 
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log(P orb [d]) 

Figure 8. Age of the binary versus orbital period for Jcmb 
(bottom panel) and Jrmb (top panel). The open squares mark 
the positions of the PCEBs immediately after the CE-phase 
(Pce, ievoi) whereas the filled squares indicate the positions 
when the mass transfer starts (P s( j, icvoi + t coo \ + i sd ). The 
solid lines show the orbital evolution of EG UMa and HR Cam. 
The dashed horizontal line indicates the age of the universe 
t ~ 1.3 x 10 10 yrs (see text). 



Fig. 6; see also Fig. 8), the number of pre-CVs among the 
PCEB sample is nearly independent on it (see Table 2, 
Fig. 8). From the sample of pre-CVs 5 (3) will initiate 
mass transfer at orbital periods > 3h whereas the other 
9 (13) systems will start their CV life in or below the or- 
bital period gap. This finding is in good agreement with 
the prediction of King et al. (1994) that >67% of CVs will 
start mass transfer near or below the period gap. Possible 
observational selection effects will be discussed in Sect. 9. 

A further test on whether the current pre-CV pop- 
ulation is a representative sample of progenitors for the 
present CV population comes from comparing their space 
density to that of CVs. An estimate of the space density 
of PCEBs/pre-CVs requires obviously the knowledge of 
their distances. 



7. Distance estimates 

We estimate the distances to the systems in Table 1 for 
which no previous distance determination has been avail- 
able in the literature using the following methods. 



7.1. PCEBs with a white dwarf primary 

For all but one system without a published distance esti- 
mate ultraviolet (UV) spectroscopy is available in the pub- 
lic archives of the International Ultraviolet Explorer (IUE) 
and the Hubble Space Telescope. We fitted the UV data 
listed by Table 3 with pure-hydrogen white dwarf model 
spectra computed with the code described by Gansicke 
et al. (1995), fixing \ogg = 8.0, i.e. M x « 0.6 M & . The 
free parameters of this fit were, hence, the white dwarf 
temperature T c ff and the scaling factor f/H — inRi 2 /d 2 , 
with / the observed flux and H the Eddington flux of 
the model spectra. The results from our fits are given in 
Table 3. Overall, the white dwarf temperatures that we 
derived here are in good agreement with the published 
values. 

Fixing the surface gravity ( = white dwarf mass) in the 
fit introduces a systematic uncertainty in the derived dis- 
tances, as the white dwarf radius obviously depends on the 
assumed mass. However, it is important to notice that also 
the white dwarf temperature derived from the fit depends 
on the assumed mass, with higher (lower) temperatures re- 
sulting for higher (lower) log<?. This effect compensates to 
some extent the R\(Mi) dependence, and the error in the 
spectroscopic parallaxes is ~ ±15% for an assumed range 
\ogg = 8.0 ±0.5 (i.e. 0.35 Mq £ M Y <, 0.9 M©), which 
covers probably most of the analysed objects. The Hyades 
member HZ 9 can be used to estimate the robustness of 
our distance estimate and error analysis. From our model 
fit to unpublished HST STIS spectroscopy, we find a tem- 
perature which is slightly below the low end of the range 
quoted by Guinan & Sion (1984), which was based on the 
analysis of IUE spectroscopy. However, the much better 
HST data available now clearly shows the 1400 A quasi- 
molecular H2+ absorption which indicates a temperature 
< 20 000 K (Fig. 9). We derive from the scaling factor of 
the model a distance of d = 40±6 pc, which is entirely con- 
sistent with the Hipparcos-measured parallax (Perryman 
et al., 1998). 

For GK Vir, we had to fix the temperature to the value 
derived by Fulbright et al. (1993) as the temperature- 
sensitive Lya absorption line is entirely blended with geo- 
coronal emission in the available IUE spectrum. 

A final note concerns the distance of LM Com. Orosz 
et al. (1999) quote a distance of 290-308 pc, which they 
base on a lengthy discussion of the white dwarf con- 
tribution to the observed R band magnitude. However, 
as the contamination of the M-dwarf is lowest in the 
blue, we re-examined the distance of LMCom scaling 
a T cff = 30 000 K, log 5 = 8.0 model to F A (4700l) ~ 
17 x 10~ 16 erg cm _2 s _1 A (measured from Orosz et al. 
1999 Fig. 7), and find d = 170 pc. 

7.2. PCEBs containing a sdB primary 

Considering the rather similar properties of the sdB 
primaries in HWVir, HS 0705+6700, PG 1017-086, and 
NY Vir (Mi, T e g, see Table 1), we computed (crude) dis- 
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Figure 9. The HST/STIS echelle spectrum of HZ9 along with 
the best-fit T e g = 17 400 K pure hydrogen white dwarf model, 
assuming logp = 8.0. The sharp "absorption dips" at the red 
end of the spectrum are due to gaps between the echelle orders. 

Table 3. Distance estimates for PCEBs, see text 



Ppceb - 2.9x10" 



0.5 



PCEB 



6.6 x 10- V 
2.9 x 10~ 5 pc 



(Jcmb) 
(Jrmb) 



AML prescriptions results from the fact that the mean 
evolution time scale is significantly longer for J = Jrmb 
and, thus, the existence of more old and yet undetected 
PCEBs is predicted. However, even the higher value of 
Ppceb, assuming Jrmb, is still significantly below the the- 
oretical predictions of de Kool & Ritter (1993, their Fig. 4) 
indicating that our sample is probably not only biased to- 
wards young systems (see Sect. 9.3). 

Requiring itot < to for PCEBs which can be con- 
sidered pre-CVs, we find 6 (Jcmb) and 7 (Jrmb) of 
these systems within 100 pc, resulting in lower limits on 
their space density yO prc cv > 1-4 x 10~ 6 pc~ 3 (Jcmb) and 
PpreCV > 1-7 x 10~ 6 pc~ 3 (Jrmb)- We estimate the actual 
pre-CV space density, again taking into account age and 
evolution time scale and find 



PprcCV — 



1.7 x 10~ 6 pc~ 3 
6.0 x 10~ 6 pc~ 3 



(Jcmb) 
(Jrmb). 



(13) 



Object 


Instrument 


Dataset 


T efI [K] 


d[pc] 


EG UMa 


HST/FOS 


y 16u0502t 


13 400 


32 


HR Cam 


HST/STIS 


o6gj02020 


20 800 


72 


UZ Sex 


IUE 


swp27393 


17 200 


35 


BPM 6502 


IUE 


swp27351 


21400 


25 


HZ 9 


HST/STIS 


o5dma6010 


17400 


40 


CC Cet 


IUE 


swp27392 


25 000 


89 


GK Vir 


IUE 


swp07459 


48 800* 


350 



Our estimate of the current pre-CV space density 
allows an estimate for the present CV space density. 
Considering the average of the pre-CV lifetime (i pre cv) of 
the systems we find within 100 pc and following Politano 
(1996) in assuming that the CV birthrate has been con- 



stant since t ga \ = 10 10 yrs we get 



Pew 



, PprcCV , 1n _5 
Igal- ~ 1 X 10 



pc 



(14) 



* Temperature fixed to the value given by Fulbright et al. 
(1993). 

tances estimates by scaling the distance of HW Vir (d = 
171 pc, V = 10.6; Wood & Saffer 1999) to the observed 

V magnitudes of the other sdB PCEBs (HS 0705+6700: 

V = 14.8; PC 1017-086: V = 14.4; NY Vir: V = 13.2). 
This estimate ignores possible extinction as well as differ- 
ences in T e ff and R\ . 



8. Space densities 

Inspection of Table 2 shows that 12 PCEBs from our sam- 
ple are within 100 pc (including RE 1016-053, d = 75 - 
109 pc). This gives an "observed" space density of (= 2.9 x 
10~ 6 pc~ 3 ) which can be considered to be a conservative 
lower limit on the PCEB space density as the presently 
known PCEB sample is systematically dominated by 
young PCEBs (see Sect. 5.1, Table2). Accounting for this 
bias, we estimate the true PCEB space density to be: 



(12 



where /pceb is the mean fractional PCEB life time that 
the systems have already passed through. The substan- 
tial difference between the estimates of ppceb for the two 



ipreCV 

for both angular momentum loss prescriptions. 

Considering the assumptions involved and the small 
number statistics of the known pre-CV sample, this re- 
sult has to be regarded as a rather rough estimate. 
Nevertheless, our estimate for pev is in agreement with 
the high end of current observational estimates (see 
Gansicke et al., 2002) but below theoretical predictions: 
10- 4 pc- 3 (dc Kool, 1992) to 2 x 10~ 5 pc~ 3 (Politano, 
1996). As we will discuss in Sect. 9.4, the presently known 
sample of PCEBs/pre-CVs is very likely incomplete not 
only with respect to those systems containing cold white 
dwarfs, but also to those systems containing an early-type 
secondary. Hence, we would like to stress that all the space 
densities derived above should be regarded as lower limits 
to the true values. 

Considering the binary age postulate (BAP) of King 
& Schcnker (2002) outlined in Sect. 2 we note that assum- 
ing reduced magnetic braking indeed leads to an increase 
of the averaged evolution time scale for pre-CVs but the 
presently known PCEB population is too small and too 
strongly biased to prove or disprove the BAP scenario. 

9. Observational selection effects 

We have shown that the known population of PCEBs is 
dominated by young hot systems (Sect. 5.1), and that, ap- 
plying the current theories of (pre-)CV evolution, the ma- 
jority of them will remain in a detached configuration for 
many 10 8 yrs to several 10 9 yrs (Sect. 6) - sufficiently long 
for their white dwarf primaries to cool - depending on 
t sd and Mi - to temperatures 6000 K & T cS <> 15 000 K 
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Figure 10. The currently known populations of single white dwarfs (top) and PCEBs (bottom) as a function of their effective 
temperature. Here, we follow McCook & Sion (1999) and plot the spectral temperature index, O of f = 50 400/T e fr, which is a 
direct indicator of the age of the white dwarfs. The grey shaded histograms represent systems which have been found in the 
PG survey. Obviously, there exists a population of cool single white dwarfs whereas PCEBs white dwarfs are generally hot and 
young. 



(see Fig. 7) These findings strongly suggest that - if the 
reality about AML in PCEBs is somewhere bracketed be- 
tween the prescriptions of Jcmb and Jrmb _ the currently 
known population of PCEBs is highly incomplete. 

As a test for possible selection effect, we compare the 
effective temperature distribution of the white dwarfs in 
our PCEB sample to that of field white dwarfs. 



9.1. The known PCEB population vs. the McCook & 
Sion (1999) white dwarf population 

McCook & Sion (1999) present a catalogue of 1793 spec- 
troscopically identified white dwarfs, including informa- 
tion on the white dwarf temperatures where available. 
The temperature distribution of the McCook & Sion 
white dwarf sample peaks at temperatures T c k vvd — 
15 000 - 22 000K (top panel of Fig. 10) with a mean ef- 
fective temperature of T c s wd = 19 360 K. Also shown 
in Fig. 10 (bottom panel) is the temperature distribu- 
tion of the PCEB white dwarf primaries from our sam- 
ple (Table 1), which also peaks in the range T e ff pceb — 
15 000 — 22 000 K, but lacks the low-temperature tail seen 
in the McCook & Sion sample. The average temperature 
of the white dwarfs in our PCEB sample is T e g pceb = 
33 423 K. Using \ 2 statistics, the significance that two dis- 
tributions arc different is 



P ( v - ^ = /o X/V 



r 2 J / °° f/2-i e -*dt 



(15) 



where v is the number of degrees of freedom, i.e. the num- 
ber of bins, and x 2 is given by: 



(16) 



with Ni the number of objects in bin i, N w( \ = X^=i -^wd,i 
and N pc = Yli=i -Wpc.i- We obtain for the probability that 
the two distributions are different: 



PI - 

1 2' 2 



0.93. 



(17) 



It appears, hence, likely that the PCEB white dwarfs 
are indeed on average hotter - and therefore younger - 
than the single field white dwarfs of (McCook & Sion, 
1999). 

9.2. PCEB and field white dwarfs in the Palomar 
Green survey 

Whereas the McCook & Sion (1999) is a heterogenous 
catalogue of white dwarfs identified by various means, it 
is of fundamental interest to compare the properties of 
white dwarfs in PCEBs and of field white dwarfs drawn 
from a single survey. Inspection of Table 1 shows that 
five of our PCEBs with a white dwarf primary are con- 
tained in the Palomar-Green (PG) survey (Green et al., 
1986): UZSex, CCCct, LM Com, GKVir and NNSer. 
Fig. 10 shows the temperature distribution of both the 
field white dwarfs (upper panel) and of the PCEB white 
dwarfs (bottom panel) from the PG survey. The mean val- 
ues of the effective temperatures for PCEB white dwarfs 
and single white dwarfs contained in the PG survey are 
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^effpcEB = 35 380 K and T cff WD ~ 21 000 K, respec- 
tively. The mean effective temperature of the white dwarfs 
in PCEBs is again somewhat higher than that of their field 
relatives. As there are only five PCEBs with a white dwarf 
primary included in the PG-survey, the statistic is very 
poor but for completeness we give the probability that 
the difference of the two distributions (see the shaded his- 
tograms in Fig. 10) can be explained by chance 

1 ~ p t) = °- 5L (18) 

Obviously, from this subset of systems we can neither say 
that the white dwarfs in PCEBs are systematically hotter 
and younger than their field relatives nor can we exclude 
that this is the case. In the next section we will discuss ob- 
servational selection effects. This will give us a hint which 
possibility is more likely. 

Notice, the probabilities given by Eq. (17) and (18) are 
calculated by neglecting the very few detected extremely 
old single white dwarfs, i.e. using only 9 (Eq. (17)) respec- 
tively 5 (Eq. (18)) bins. 

9.3. Finding PCEBs 

Based on the presently available data (Sect. 9.1 and 9.2), 
we can not exclude that the white dwarfs in PCEBs are 
- on average - systematically hotter (and younger) than 
single white dwarfs. This statement may be entirely based 
on selection effects skewing the temperature/age distribu- 
tion of the white dwarfs in the known PCEB population. It 
is, however, very difficult to assess this bias quantitatively. 
We list in Appendix A for the PCEBs from our sample in 
Table 1 (1) the way that they were discovered in the first 
place, and (2) the way that their binarity has been identi- 
fied. Inspection of these notes shows that the vast majority 
of the known PCEBs have been identified in the first place 
because of their white dwarf component - either as blue 
object or X-ray source. Obviously, the presence of a late- 
type companion may significantly affect the colour of a 
PCEB: the cooler the white dwarf and/or the earlier the 
companion are, the less blue the system will appear. This 
may affect different candidate samples in different ways, 
depending on how the criterion blue has been defined. 

We have computed the colours expected for PCEBs 
containing a white dwarf with 6000 K< T cff < 50 000 K 
and an ZAMS companion with spectral type K0 to M6, 
the results for U — B and B — V arc shown in Fig. 11. 
Comparing these simulated colours to the selection cri- 
terion of the PC survey, U — B < 0.46 (Green et al., 
1986), it becomes apparent that PCEBs containing a cool 
(T cff <, 15 000 K) white dwarf will be included in the PG 
survey only if their companions are of spectral type ~M4 
or later. Similarly, for such PCEBs B — V > 0.5, which 
will most likely prevent their selection as white dwarf can- 
didates. As an example for the "historic" white dwarf 
searches, Giclas et al. (1965) used a Lowell colour class 
-1 or (corresponding to B - V < 0.0 or B - V < +0.2, 



and U - B < -0.78 or U - B < -0.60) as selection cri- 
terion. Consequently, the GD lists will contain only those 
PCEBs containing moderately hot white dwarfs and late 
type secondaries, such as HRCam and MS Peg. 

PCEBs which were not initially selected as white dwarf 
candidates are: (1) the planetary nebulae; (2) high proper 
motion objects (RRCae, BPM6502, and BPM 71214) - 
interestingly enough, RR Cae is the PCEB containing the 
coldest white dwarf in our sample; and (3) variable stars 
(V471Tau: spectroscopic; BEUMa: photometric). In the 
case of V471 Tau the optical emission is largely dominated 
by the K2V companion, and with B~V ~ 0.76 the system 
does clearly not qualify as a white dwarf candidate. 

9.4. The lack of early companions 

So far we have focussed much of our attention on the tem- 
peratures of the white dwarfs in our PCEB sample, as 
this parameter is a direct measure of the age of these sys- 
tems. However, Table 1 suggests that this sample is also 
biased with respect to the mass of the companion stars. 
The average mass of the companion in our PCEB sample 
is M scc pceb = 0.24, whereas CVcat (Kube et al., 2002) 
lists the masses for 80 CVs with M sec cv = 0.43. A di- 
rect consequence of this dominance of low mass companion 
stars in the known PCEB sample is that the majority of 
these systems evolve into contact at short orbital periods 
(Sect. 6; Fig. 5). 

The shortage of PCEBs with massive/early- type donor 
stars (0.4M Q <> M sec <> 0.8 M ), i.e. the progenitors of 
long-period CVs, is subject to the selection effect which 
we just discussed in Sect. 9.3: in PCEBs with an early-type 
companion the secondary contributes significantly to the 
total optical emission of the system, and, consequently, 
such binaries do not qualify as white dwarf candidates 
(Fig. 11), the major discovery channel for PCEBs. 

As a result of the (very likely) incompleteness of known 
PCEBs with an early-type donor the space density esti- 
mates presented in Sect. 8 should be considered to be lower 
limits to the true ppceb and p pre cv- 

10. Conclusion 

We have presented the properties of the currently known 
sample of well-observed PCEBs and have calculated their 
past and future evolution based on two different assump- 
tions for the angular momentum loss mechanism, adopting 
either the "classical" disrupted magnetic braking prescrip- 
tion (Jcmb) or the AML prescription derived by Sills et al. 
(2000) (Jrmb)- The results from this study are: 

1. The presently known PCEB population is dominated 
by young systems most of which have evolved only 
through a small fraction of their lifetime as detached 
binaries. On average, the white dwarfs in PCEBs are 
hotter than field white dwarfs. 

2. While the evolution of an individual system strongly 
depends on whether we assume JRMBOr Jcmb, the 
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Figure 11. Simulated colours for PCEBs containing a white dwarf with 6000 K < T e g < 50 000 and a companion with a 
spectral type in the range K0 to M6. The left panel shows in gray the PCEB parameter space that do not fulfill the selection 
criterion of the PG survey (U — B < 0.46). 



total number of pre-CVs within the PCEB sample, 
i.e. systems which will start mass transfer within a 
Hubble time, is nearly the same for both prescriptions. 
Considering the time scale on which the PCEBs will 
evolve into semi-detached CVs we predict the existence 
of a large population of old PCEBs containing cold 
white dwarfs. The present lack of such systems is very 
likely the result of observational biases, as the major- 
ity of the known PCEBs were initially selected as blue 
objects. 

3. An additional consequence of the observational selec- 
tion effects involved in the discovery of PCEBs is a 
shortage of systems containing a "massive" /early-type 
companion (0.4 M ^ M scc 0.8 M ) - which are 
the progenitors of long-period CVs. Indeed, our calcu- 
lations predict that most of the pre-CVs among the 
PCEBs in our sample will evolve into CV with short 
orbital periods (P or b < 3 h) . 

4. The space density of PCEBs estimated from the cur- 
rently known sample is 6 x 10~ 6 pc~ 3 <ppcEB~3 x 
10 _5 pc -3 , depending on the assumed AML prescrip- 
tion. Taking into account the age as well as the evo- 
lution time scale of the pre-CVs in our sample and 
assuming that the birthrate of PCEBs remains un- 
changed we derive for the CV space density pcv ~ 
10 _5 pc~ 3 , nearly independent on whether we assume 
Jcmbot Jrmb- This value is somewhat higher than 
current observational estimates but below theoretical 
predictions. However, due to the observational bias 
working against the discovery of PCEBs containing an 
early-type companion our estimates should be consid- 
ered to be lower limits to the true values of pcv and 
Ppceb- 
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Note added in proof: After the submission of our paper, 
Raymond et al. (2003) have published a list of 109 PCEBs 
found in the Sloan Digital Sky Survey. The average white 
dwarf temperature of this sample 16 000 K, significantly 
lower than that of our sample of PCEBs with known or- 
bital period (Table 1), and demonstrates that more sophis- 
ticated selection criteria (compared to pure "blue" sur- 
veys) can successfully identify PCEBs containing rather 
cold white dwarfs. Unfortunately, only a single system in 
the list of Raymond et al. (2003) has an orbital period 
measurement - determining the periods of their PCEBs 
harbouring white dwarfs with T e g ^ 12000 K provides a 
direct test for our prediction (2) above. 

Appendix A: The discovery history of our PCEB 
sample (Table 1) 

RR Cae: Discovered as high proper motion object by 
Luyten (1955), DA white dwarf plus Balmer emission 
line spectrum noted by Rodgers & Eggen (1974), late- 
type nature for the companion suggested by Bessell & 
Wickramasinghc (1979), discovery of eclipses and orbital 
period measurement (photometric) by Krzeminski (1984). 
Two comprehensive studies by Bruch & Diaz (1998) and 
Bruch (1999). 

EG UMa: Discovered as a white dwarf on Schmidt 
prism plates by Stephenson (1960), emission lines noted 
by Greenstein (1965), orbital period (spectroscopic) by 
Lanning (1982). Recent comprehensive study by Bleach 
et al. (2000). 

EC 13471-1258: Discovered in the Edinburgh-Cape 
faint blue object survey of high galactic latitudes, bina- 
rity revealed through the detection of eclipses (Kilkenny 
et al., 1997). Recent study by Kawka et al. (2002). 

BPM 71214: Discovered as a high proper motion object 
(Luyten, 1963). Recent study by Kawka et al. (2002). 
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HR Cam: Listed as WD candidate by Giclas et al. 
(1970), spectroscopic confirmation of WD nature by 
Wills & Wills (1974), detection of a red companion by 
Zuckerman & Becklin (1992), orbital period (spectro- 
scopic and photometric) from the comprehensive studies 
of Marsh & Duck (1996); Maxted et al. (1998). 

UZSex: Discovered in the PG survey, listed as a 
DA/composite spectral type by Green et al. (1986), or- 
bital period (spectroscopic) from the comprehensive study 
of Saffer et al. (1993), see also the recent analyses by Bruch 
& Diaz (1999); Bleach et al. (2000). 

BPM 6502: Discovered as a high proper motion star by 
Luyten (1957a), listed as WD candidate by Eggen (1969), 
spectroscopically confirmed as a WD by Wegner (1973), 
detection of companion and orbital period (spectroscopic) 
by Kawka et al. (2000). 

HZ 9: Discovered in a search for faint blue stars by 
Humason & Zwicky (1947), mentioned as a DA+dMe 
candidate by Greenstein (1958), orbital period (spectro- 
scopic) measured by Lanning & Pesch (1981). 

MS Peg: Listed as WD candidate by Giclas et al. 
(1965), spectroscopic confirmation of WD nature by 
Greenstein (1969), emission lines and radial velocity varia- 
tions detected by Tytler & Rubenstein (1989) and Schultz 
et al. (1993), orbital period (spectroscopic) measured by 
Schmidt et al. (1995). 

HS 0705+6700: Classified in the Hamburg Quasar 
Survey (Hagen et al., 1995) as hot star candidate, eclipses 
discovered and orbital period measured by Drechsel et al. 
(2001). 

HWVir: Listed as ultraviolet-bright star by 
Carnochan & Wilson (1983) and classified as likely 
sdB star by Berger & Fringant (1980), discovery of 
eclipses and orbital period (photometric) by Menzies 
(1986). 

LM Com: Listed as blue object by Iriarte & Chavira 
(1957), also detected in the PG survey (Green et al., 1986), 
and identified as DA-fdM binary by Ferguson et al. (1984), 
Ha emission detected by Orosz et al. (1997), orbital period 
(spectroscopic) and comprehensive study by Orosz et al. 
(1999). 

PG 1017-086: Discovered in the PG survey, classified 
as sdB star by Green et al. (1986), variability discovered 
and orbital period measured by Maxted et al. (2002). 

CC Get: Discovered in the PG survey, listed as a 
DA/composite spectral type by Green et al. (1986), or- 
bital period (spectroscopic) from the comprehensive study 
of Saffer et al. (1993), see also Somers et al. (1996). 

V471 Tau: The brightest and best-studied pre-CV. 
Listed as a spectroscopic binary by Wilson (1953), discov- 
ery of eclipses, spectral classification as DA+dK binary, 
and measurement of the orbital period (photometric) by 
Nelson & Young (1970). 

PG 1336-018: Discovered in the PG survey, classified 
as sdB star by Green et al. (1986), eclipses discovered and 
orbital period measured by Kilkenny et al. (1998). 

AA Dor: Listed as faint blue star by Luyten (1957b) 
and as spectroscopically variable object in the foreground 



of the LMC by Feast et al. (1960), discovery of eclipses 
and measurement of the orbital period (photometric) by 
Kilkenny et al. (1978). 

RE2013+400: Discovered as bright EUV source dur- 
ing the ROSAT Wide Field Camera all-sky survey and 
listed as a WD by Pounds et al. (1993), variable Balmcr 
emission was noted by Barstow et al. (1993), the orbital 
period (spectroscopic) was measured by Thorstensen et al. 
(1994). 

GK Vir: Discovered in the PG survey, eclipse discov- 
ered and orbital period (photometric) measured by Green 
et al. (1978). 

MTSer: Classified as PN on the Palomar Sky Survey 
plates by Abell (1955), variability suggested by Abcll 
(1966), recurrent photometric variability was detected by 
Grauer & Bond (1983) who also suggested the sdO+dM 
nature of the system, for a more comprehensive photomet- 
ric study, see Bruch et al. (2001). 

IM CMa: Discovered as bright EUV source during the 
ROSAT Wide Field Camera all-sky survey and listed as 
a WD by Pounds et al. (1993), Balmer emission detected 
and orbital period (spectroscopic) measured by Vennes & 
Thorstensen (1994). 

NNSer: Discovered as a CV candidate in the PG sur- 
vey by Green et al. (1986), discovery of deep eclipses, clas- 
sification as DA+dM binary and measurement of the or- 
bital period (photometric) by Hacfncr (1989). 

RE 1016-053: Discovered as bright EUV source during 
the ROSAT Wide Field Camera all-sky survey and listed 
as a WD by Pounds et al. (1993), classified as DA+dM 
by Jomaron et al. (1993), orbital period (spectroscopic) 
measured by Tweedy et al. (1993). 

V477 Lyr: Classified as PN on the Palomar Sky Survey 
plates by Abell (1955), variability suggested by Abell 
(1966), eclipse discovered and period (photometric) mea- 
sured by Bond (1980). 

Abell 65: Classified as PN on the Palomar Sky Survey 
plates by Abell (1966), an estimate of the orbital period 
(photometric) was published by Bond & Livio (1990). 

KVVel: Listed as luminous star in the list of 
Stephenson & Sanduleak (1971), planetary nebula na- 
ture suggested by Holmberg et al. (1978), spectroscopi- 
cally identified as sdO star within a planetary nebula by 
Drilling (1983). Photometric variability discovered and or- 
bital period (photometric) measured by Drilling & Bravo 
(1984). 

UU Sge: Classified as PN on the Palomar Sky Survey 
plates by Abell (1955), variability suggested by Abcll 
(1966), eclipses discovered and orbital period (photomet- 
ric) measured by Miller et al. (1976). 

HS 1136+6646: Classified in the Hamburg Quasar 
Survey (Hagen et al., 1995) as hot star candidate, iden- 
tified as a spectroscopic binary by Heber et al. (1996), 
orbital period (spectroscopic) measured by Sing et al. 
(2001). 

BE UMa: Identified as variable star and orbital period 
(photometric) by Kurochkin (1964, 1971), recovered as 
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emission line star in the PG survey by Ferguson ct al. 
(1981). 
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